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Abstract 

We describe an algebraic algorithm which allows to express every one-loop 
lattice integral with gluon or Wilson-fermion propagators in terms of a small 
number of basic constants which can be computed with arbitrary high precision. 
Although the presentation is restricted to four dimensions the technique can 
be generalized to every space dimension. Various examples are given, including 
the one-loop self-energies of the quarks and gluons and the renormalization 
constants for some dimension-three and dimension-four lattice operators. We 
also give a method to express the lattice free propagator for Wilson fermions 
in coordinate space as a linear function of its values in eight points near the 
origin. This is an essential step in order to apply the recent methods of Liischer 
and Weisz to higher-loop integrals with fermions. 

PACS: 11.15.Ha 12.38.Gc 



1 Introduction 



Perturbation theory plays an important role in our present understanding of quan- 
tum field theory. In particular on the lattice Feynman diagram computations are 
performed to obtain such quantities as ratios of A-parameters, non-universal coeffi- 
cients of /3-functions or renormalization constants of lattice operators. Due to the 
loss of Lorentz invariance lattice calculations are usually particularly involved and 
thus in order to obtain reliable results one has to make use of computer symbolic 
programs. To implement this strategy an important step is the simplification of the 
lattice integrals appearing in the Feynman diagrams of the theory. In we presented 
a general technique which allows to express every one-loop bosonic integral at zero 
external momentum in terms of two unknown basic quantities which could be com- 
puted numerically with high precision. This method allows the complete evaluation 
of every diagram with gluon propagators. 

In this paper we want to generalize the technique to deal with integrals with both 
gluonic and fermionic propagators. For the fermions we use the Wilson action [Q. 
Notice that our method depends only on the structure of the propagator and thus it 
can be applied in calculations with the standard Wilson action as well as with the 
improved clover action 0]. We show that every integral at zero external momentum 
can be expressed in terms of a small number of basic quantities (nine for purely 
fermionic integrals, thirteen for integrals with bosonic and fermionic propagators). 
The advantage of this procedure is twofold: first of all every Feynman diagram can be 
computed in a completely symbolic way making it easier to perform checks and verify 
cancellations; moreover the basic constants can be easily computed with high precision 
and thus the numerical error on the final result can be reduced at will. Although the 
presentation is restricted to four dimensions the technique can be generalized to every 
space dimension. 

A second important application of our method is connected with the use of 
coordinate-space methods for the evaluation of higher-loop Feynman diagrams. This 
technique, introduced by Liischer and Weisz @], is extremely powerful and allows a 
very precise determination of two- and higher-loop integrals. One of the basic ingredi- 
ents of this method is the computation of the free propagator in coordinate space. We 
will give an algorithm which allows the analytic determination of the free fermionic 
propagator in x-space in terms of eight basic constants which can be reinterpreted as 
the values of the propagator near the origin. 

The paper is organized as follows: in section 2 we review the computation of the 
continuum limit of lattice integrals, showing that the calculation can be split in two 
parts: the evaluation of a subtracted continuum integral and of a certain number 
of lattice integrals with zero external momentum which are then discussed in the 
following sections. In section 3 we present the method in the bosonic case, simplifying 
the strategy discussed in fTj]. Section 4 represents the core of the paper and gives 
the detailed algorithm for fermionic and mixed bosonic-fermionic integrals. In both 
these sections we need to introduce an infrared cut-off to regularize the integrals for 
k = 0. We choose to introduce a mass m; other possibilities are discussed in section 5 
where the connection with dimensionally regularized integrals is presented. Finally in 
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section 6 we present a few examples: we give the analytic expressions for the lattice 
gluon and fermion self-energy, for the renormalization constants of dimension-three 
bilinear fermion operators and a computation of the renormalization constants for 
the operators which show up in the energy-momentum tensor and which are relevant 
in the computation of the structure functions which appear in the deep-inelastic 
scattering |0. Finally we discuss an algorithm for the computation of the free lattice 
propagators. 



2 Continuum limit of lattice integrals 

In this Section we want to discuss the computation of one-loop Feynman diagrams 
on the lattice. In general each graph has the form 

G{{p.}) = J (2.1) 

where {pi} is a set of external momenta. Of course one is not interested in the exact 
computation of but only in its value in the continuum limit. If the integral is 
ultraviolet-convergent one can simply substitute {pi}) with its continuum coun- 
terpart and obtain a continuum convergent integral. Let us now suppose that (|2J] ) 
is divergent and, for simplicity, that there is only one external momentum. If every 
propagator is massive so that F{q; {pi}) is finite for any set of momenta going to zero, 
one can use the general technique of BPHZ |^ which have been generalized on the 
lattice by Reisz writing 

G{p) = / [F{q;p) - {T--F){q;p)] + J ^.(T^- F){q; p) 

^ G%p)+G'^{p) (2.2) 
where np is the degree of the divergence of the integral and 



(T"-F)(g;p) = ^-p,,...p,, 
fc=o ■ 



9 d ' 

dpui op. 



(2.3) 

^Ml ^J^/J-k J p=o 

The first integral in (|2.2| ) is ultraviolet-finite and thus one can take the continuum 
limit. Thus all the effects of the lattice regularization remain only in the second term, 
which is simply a polynomial in the external momentum with coefficients given by 
lattice zero-momentum integrals. 

If the integrand contains massless propagators one has to be more careful: indeed 
an expansion around p = can give rise to infrared divergences. A simple way 
out consists in introducing an intermediate infrared regularization: one can use the 
dimensional regularization P, |l]] working in dimension c? > 4, or introduce a mass in 
the propagators. In both cases G^{p) and G^{p) will be singular for d ^ 4 or m — 
but, of course, the singularity will cancel when summing up the two terms. 

In conclusion the computation of the continuum limit of ( |2.1D splits in the compu- 
tation of two different quantities: a continuum ultraviolet-finite integral and a certain 
number of zero-momentum lattice integrals whose computation will be discussed in 
the next sections. 
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3 Bosonic integrals 

In this section we discuss the evaluation of the most general one-loop lattice integral 
at zero external momentum with bosonic propagators. It is very easy to see that any 
such integral can be written as a linear combination of terms of the form 

B{p- n., n„ n., n,) = j ^ ^ \l^k^'my ^^'^^ 

where p and rii are positive integers, = 2sin(fc^/2) and 

DB{k,m) = k'^ + . (3.2) 

is the inverse bosonic propagator. In the following when one of the arguments is 
zero it will be omitted as argument of B. 

We will review here a general technique for expressing these integrals in terms of 
three constants 

We will first generalize ( |3.1| ) by considering the following more general integrals 

, r d^k fcf 4f 4f 4f* ^ ^ 

Bsip; n.., n„ n„ n,) = ^ Ds{k,mY^^ ^'-'^ 

where p is an arbitrary integer (not necessarily positive) and 5 a real number which is 
introduced in order to avoid singular cases at intermediate stages of the computation 
and which will be set to zero at the end. 

The first result we want to prove is that each integral Bs{p]nsc,ny,nz,nt) can be 
reduced through purely algebraic manipulations to a sum of integrals of the same 
type with rix = riy = Uz = rit = 0. 

Indeed the integrals Bs satisfy the following recursion relations: 

Bs{p;l) = ^[Bs{p-l)-m'Bs{p)] 

Bs{p;x,l) = ^[Bs{p-l;x) - Bs{p]x + 1) -m'^Bs{p;x)] 

J3s{p;x,y,l) = ^[Bs{p ~ l;x,y) - Bsip;x + l,y) - Bs{p;x,y + 1) 

~m^Bs{p] x,y)] 

Bsip;x,y,z,l) = Bsip - I; x,y, z) - Bsip; x + l,y, z) - Bsip; x,y + 1, z) 

-Bsip; a;, ^ + 1) - m^Bsip; ^) (3-4) 
which can be obtained using the trivial identity 

4 

DBik,m) =Y,kj+m^ . (3.5) 



Furthermore, when r > 1, we can write 



iklr , Ckly-' ^2-SplsinA;^-^— — -i— — (3.6) 



DBik,m)P+^ DBik,m)P+^ p + 6-1 dk^ DBik,m)P+^-^ 
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Then, integrating by parts, we obtain the recursion relation: 



Bs{p] . . . ,r) 



r — 1 4r — 6 

-Bs{p-l;...,r-l) r^—Bs - 1; . . . , r - 2) 



p + ' ' ' p + 6-1 

+ 4Bs{p;...,r-l) (3.7) 

Let us notice that for p ^ 1 this recursion relation is regular for 5 — » 0. For p = 1 
instead the coefficients of Bs{0; . . .) diverge as 1/6. This means that to compute 
Bs{l; . . .) for (5 = 0, we need to compute Bs{0; . . .) including terms of order 6. Since by 
the application of the previous recursion relations to Bs{0] . . .) we generate Bs{—1] ■ ■ ■), 
then Bs{—2] . . .) and so on and their coefficients are finite for 5 — >■ 0, we see that in 
general we need to compute all integrals Bs{p; n^, Uy, n^, rit) with p < up to 0(5^). 

The previous relations allow to reduce every integral Bsip; rix, Uy, Uz, rit) to a sum 
of the form 

p 

Bs{p]n^,ny,nz,nt) = ^ ar{'m,6)Bs{r) (3.8) 

r=p-nx-ny-nz—nt 

The m-dependence of ar{m,6) is very simple: it is indeed a polynomial in m^. Let 
us now discuss the 5-dependence. If p < only nonpositive values of r are allowed in 
( p.8|) with coefficients which are regular for 6 0. For p > the situation is more 
complicated: for r > 1, lim^^o ar(^, 6) is finite while for r < ar{m, 6) may behave 
as 1/(5 when 6 goes to zero. As we already observed, this means that we need to 
compute Bs{r), r < 0, including terms of order 6. 

Now let us show that all Bs{p) can be expressed in terms of a finite number of 
them. Although this can be shown for generic values of the massQ many simplifications 
occur if one restricts the attention to the massless case, i.e. if one considers the limit 

— > keeping only the non-vanishing terms. Let us thus discuss the limit — > 
ofBsip). 

Simple power-counting shows that the integrals Bs{p), p < 1, are finite. To com- 
pute the divergent part of Bs{p), p > 2, let us start from the well-known representation 

BAp) = —-4 ;r / dW^+^-^e-"^'^'^-^^ UXf (3.9) 

2P+^T{p + 5) Jo ^ ' ^ ' 

where /o(A) is a modified Bessel function [IC]. We introduce constants bi that are 
defined by the asymptotic expansion for large x of /o( 

°° h- 

^o(a:)^-e-E-^ • (3-10) 

1=0 

The constants bi are rational numbers multiplied by l/vr^. Then the divergent part 
of Bs{p) is given by 

1 ^ k^2T{p + S -i) 



r(p + 5)t^ 2*(m2)p+'5-i 



^The analogous procedure in two dimensions is presented in Appendix A. 3 of 
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Since here p is positive we can set 6 = 0. However in the next section we will be 
interested also in the divergent part proportional to S. Thus expanding in S and 
neglecting finite terms, up to 0{6^), we getf\ 



'p~2 



1 hi.2^{p - i) 

r(p) ^ 2^(m2)p-^ 2PV{p) 



logm^ 



-5 



1 t^h_,T{p-i) 



V{p) 



1=2 



logm^+ ^ - 



■Jp-2 



1. 2 "^'1 



k=p—i 
2 



- log + ^ - ) log m 



k ^ 

+ 0(5^) 



2PT{p) 

Let us now go back to ( |3.8|) rewriting it as 

B5{p;n^,ny,n^,nt) = J^^riO, S)Bs{r) +J2Mfn, 5) - ar{0, 5))Bsir) 



(3.12) 



(3.13) 



and let us consider the limit —>■ 0. It is clear that in the second sum only Bs{r) with 
r > 3 can contribute, since only these integrals have power divergences for — 0. 
As the values of r in the sums satisfy r > p, we find that the second sum contributes 
only for p > 3. Thus we get for p < 2 



Bs{p;n^,ny,nz,nt) = ^ a,.(0, (5)S5(r) + 0(m^ 



(3.14) 



while, for p > 2, 



Bsip;n^,ny,n;,,nt) = '^ar{0,6)Bsir) + n{m,6)+ 0{ 



m 



(3.15) 



where TZ{m;6) is a polynomial in l/m^ whose coefficients, for 5 — 0, are rational 
numbers multiplied by l/vr^. 

Let us finally find the last recursion relations. Let us start from the trivial identity 



Bs{p; 1, 1, 1, 1) - ABsip + 1; 2, 1, 1, 1) - Bs{p + 1; 1, 1, 1, 1) 







(3.16) 



and let us apply the previous procedure to reduce each term to the form ( p. 141) and 
( ^.15| )p|. We thus get a non trivial relation involving Bs of the form 



J2 br{p;6)Bs{r)+ S{p;m,6) = 

r=p—A 



(3.17) 



^Notice that we are dealing here with two different Hmits, and (5 — + 0. In general they 

do not commute and thus it is necessary to specify the correct order in which they are taken. Here 
we first consider S —^ at fixed m and then we let m go to zero. 

■^Since we need to compute Bs{r), r < 0, including terms of order S, in the computation of the 
various terms, one must keep the contributions of order (5 if p < 0, while for p > it is enough to 
expand the identity to order 0{S). 
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where S{p; m,S) = for p < 2 while for p > 2 is a polynomial in which is finite 
for (5 — >• 0. We will use this identity to express all Bs{p) in terms of Bs{r), < r < 3. 
Indeed we can solve ( p.l7| ) in terms of Bs{p) and thus we get a relation which expresses 
it in terms of Bs{p — 1), . . . , Bs{p — 4:). We will use this relation for p > 4. On the 
other hand we can solve ( p.l7| ) in terms of Bs{p — 4:) and then shift p ^ p + A. In this 
way we obtain a relation which expresses Bs{p) in terms of Bs{p + I) , ■ ■ ■ , Bs{p + 4) . 
We use this recursion for p < — 1. Applying recursively these two relations we get 
finally (p 7^0, 1,2, 3) 

Mp) = E S)Bs{r) + r{p- m, 5) (3.18) 

r=0 

where T{p;m,6) is a polynomial in 1/m?. A direct analysis of (|3.17| ) shows the 
following properties: 

1. ifp>4, co{p;6) = 0{5% 

2. if p < -1, Cr{p] 5) = 0(5) for 1 < r < 3; 

3. if p < —1, T{p', 1TL, 6) is of order 6, while for p > 4 it is finite for 6 0. 
Substituting back in ( p.l4| ) or ( p.l5|) we get 

Bsip; n., uy, n„ rit) = A{6)BsiO) + B{5)Bsil) + C{5)Bsi2) + D{5)Bsi3) + E{m, 5) 

(3.19) 

where E{m,6) is a polynomial in 1/m^. 

We can now go back to the original integral ( |3.1| ) (notice that we are only interested 
in the case p > 0). Because of the second property of the coefficients Cr{p; S) and the 
property of T{p;m,6), we immediately see that B{6), C{6), D{6) and E{m,6) are 
finite for 5 — >■ 0. Then, as the l.h.s. is obviously finite for 5 — 0, also A{6) is finite 
for (5 — 0. Since B{0) = 1, we have finally 

B{p; n^, Uy, n„ Ut) = A{0) + B{0)B{1) + C{0)B{2) + D{0)B{3) + E{m, 0) (3.20) 

We want now to make contact with our previous work where all results were expressed 
in terms of three constants, Zq, Zi and Fq defined by 



Zi = -S(l;l,l)L^o (3.22) 



Zo = i3(l)L=o (3-21 
1 

4 

Fo = lim[167r^;B(2) + logm^ + 7£;] (3.23) 

m— >0 

It is clear how to rewrite B{1) and B{2) in terms of Fq and Zq. For B{3) a short 
calculation gives 
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Zo 0.154933390231060214084837208 
Zi 0.107781313539874001343391550 
Fn 4.369225233874758 



Table 1: Numerical values of the three constants Zq, and Fq. 



An additional simplification occurs if the original integral is finite. In this case the 
logm^ terms must cancel. They appear only in i3(2) and i3(3) and always in the 
combination (logm^ + 7^ — Fq). Thus the cancellation of logm^ implies also the 
cancellation of 7^; and Fq. All finite integrals are thus functions of and Z\ only. 
Numerical values of the constants are reported in Table |l]. 

It is interesting to notice that the same technique can be used for bosonic integrals 
in d dimensions. The basic recursions can be trivially generalized as well as the 
identity ( ^.161 ). For generic d we find that all integrals can finally be expressed in 
terms of ^^(l), . . ., Bi,{d — 1), i.e. in terms of [d — 1) constants, reducing to [d — 2) 
for infrared finite integrals. 



4 Integrals with bosonic and Wilson-fermion prop- 
agators 

We want now to discuss the computation of general integrals with fermionic and 
bosonic propagators at zero external momenta. Define 

J^{p,q;n^,ny,n^,nt) = / ^ ".^ (4.1) 



where p, q and rii are positive integers, Dsik, rrib) is defined in (|3^ ) and 



^2 



DF{k,mf) = Y,s\Yi^h + ^{k''f + m) (4.2 



is the denominator appearing in the propagator for Wilson fermionsQ. In the following 
when one of the arguments rij is zero it will be omitted as an argument of JF. Following 

^To be precise, Dp{k,mf) is the denominator in the propagator for Wilson fermions only for 
ruf — 0. For m/ 7^ the correct denominator would be 

L»F(fc,TO/) = ^sin2/ci+ (^P + m^y . (4.3) 

i 

However in our discussion nif will only play the role of an infrared regulator and thus it does not 
need to be the true fermion mass. The definition (4.2) is easier to handle than (4.3). For a discussion 
of integrals using (4.3) see Section 5.2. 
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the strategy we have used in the purely bosonic case we will first generalize the 
integrals introducing 



'Z-l 




{2n)^ DF{k,mf )P+^DB{k,mby 



(4.4) 



Here we consider both p and q as arbitrary integers (not necessarily positive). The 
parameter 5 is used in the intermediate steps of the calculation and will be set to zero 
at the end. 

To simplify the discussion we will only consider the case ryy = 1 but the technique 
can be applied to every value of rw- Moreover we will restrict our attention to the 
massless case, i.e. we will consider the integrals J^s in the hmit nib = ruf = m ^ 0. 

In the following we will present a procedure that allows to compute iteratively a 
generic J-s in terms of a finite number of them: precisely every ^^^(p, q; n^j Uy, Uz, rit) 
with q < can be expressed in terms of ^^^(1,0), ^^^(1,— 1), JF5(1,— 2), J^s{2,0), 



j^5(2,-l), ^5(2,-2), ^5(3,-2), ^5(3,-3) and ^^(3,-4); the integral ^^(2,0) ap- 



pears only in infrared-divergent integrals. If g > the result contains three additional 
constants together with the bosonic quantities Zq, Zi and Fq — 7^. 
Our procedure works in four steps: 

1. First step: we express each integral J-'s{p,q;nx,ny,nz,nt) in terms of J-'s{p,q) 



2. Second step: we express every J-'s{p,q) in terms of J-'s{r,s) with < r < 3, 
arbitrary s or r < — 1 and s = 1,2,3 or r > 4 and s — 0,-1,-2. This is 
obtained by a systematic use of the identity 



Xi(p, q) = Tsip, q- 1, 1, 1, 1) -4J^5(p, q+l- 2, 1, 1, 1) -Tn'Tsip, g+ 1; 1, 1, 1, 1) = 



3. Third step: we express the remaining Tsip.q) in terms of Ts{t,s) with r = 3, 
-4 < s < or r = 2, -4 < s < 2 or r = 1, -4 < s < 4 or r = 0, -4 < s < 6 
or r = —1 and s = 2. This is obtained by a systematic use of the identity 



X2(p, q) = Ts{p, q; 1, 1, 1, 1) - :Fsip + 1, g - 1; 1, 1, 1, 1) + J'sip + 1, q; 3, 1, 1, 1) - 



4. Fourth step: the identities Ii{p,q) and l2{p,q) are used to provide additional 
relations between the remaining integrals. We end up with the result we have 
quoted above. 

We want to notice here two basic facts. First of all, as in the bosonic case, the 
structure of the recursion relations, will force us to compute all J^sip, q] ^x, ^y, n^, rit) 
including terms of order 5 when p < 0. As a consequence we will require Xi(p, q) and 



only. 



(4.5) 



- [Ts{p + 1, ? - 2; 1, 1, 1, 1) - 2m^rs{p + 1, g - 1; 1, 1, 1, 1) 

-Fm^J^5(p+l,g; 1,1, 1,1)1 =0 



(4.6) 



9 



X2(p, q) to be satisfied up to 0(5^) respectively for p < and p < —1 and up to 0{6) 
in the opposite case, i.e. for p >1 and p > respectively. 

Let us finally mention some general properties of all the recursion relations we 
will consider: in all cases we will give results of the form 

^s{p, q;...)=Yl oipq;rs{S)J^s{r, s) + Vpg{m, 6) + O(m^) (4.7) 

rs 

where «py;rs(5) and Vpq{m, 5) will always have the following properties: 

1. if g < 0, then Q;pq;rs(5) = for s > 0. In other words we will express purely 
fermionic integrals in terms of integrals of the same type; 

2. for (5 ^ 0, we have apg-rsi^) ~ 0(1/6) for p > and r < 0, apq;rs(5) ~ 0(1) for 
p > 0, r > and p < 0, r < and apq-rsi^) ~ 0(5) for p < and r > 0; 

3. for (5 -> we have Vpq{m, S) = V^]^ (m) +0(6) for p > and q<0, Vpq{m, 5) = 
^'^pq (^) + ^i^'^) for p < and g < 0; for p > and g > we have 

Vpq{m,d) = ^{l-Slogm')V^l\m)+V^l\m) + Oid) (4.8) 

while for p < and g > we have 

Vpqim, 5) = {1- dlogm')Vi,l\m) + 6Vj,l\m) + 0{d') ; (4.9) 

in all cases V'^^{m) and V^^{m) are polynomials in 1/m^ whose coefficients are 
rational numbers multiplied by I/tt^. 

4.1 First step: the basic identities 

We will now show that in a purely algebraic way all integrals can be reduced to a sum 
of J^s{p:, q). Indeed it is easy to see that these integrals satisfy the following recursion 
relations (in giving these recursions we keep rrib ^ ruf and rw generic) 

^s{p, 1) = ^ [J^s{p, g - 1) - mlJ^sip, q) 

Ts{p,q;x,l) = ^[j^s{p,q-'i-;x) -mlJ^s{p,q;x) - Ts{p,q;x + 1) 

^s{p,q;x,y,l) = ^[j^5{p,q - l;x,y) - mlJ^s{p,q;x,y) 

-^s{p, q;x+l,y)- J^s{p, q;x,y + 1)] 

^s(p, q; X, y, z, 1) = J^s{p, q-^'^x, y, z) - mlJ^s(p, q; x, y, z) - J^s(p, q;x + l, y, z) 

-J^sip, q;x,y + 1, z) - J^s(p, q; x,y,z + 1) (4.10) 

which can be obtained from the trivial identity 

DB{k,m,) = ^^ + ml (4.11) 

i 
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A second recursion relation is obtained from the identity 

J2 H = KDb{K rub) - Dpik, ruf) - + m)) + r^(L>B(A;, m^) - mlf (4.12) 

i 

In this way we get 

•^5(p,?;2) = T5{p,q-l)-T5{p-l,q) + {m)-ml)r5{p,q) 



+^ [Fsip, g - 2) - 2mlTs{p. q-l)+ mtJ'sip, q) 



Tsip,q;x, 2) 



J^s{p, q-l;x)- Tsip - 1, q; x) 

+{m} - ml)J^s{p, q; x) - ^J^s{p, q;x + 2) 



I 'w 
3 



^siP, q-2;x)- 2mlJ^siP, q-l;x)+ mlJ^siP, q; x) 



^5{p,q;x,y,2) 



^sip, q-'i-;x,y)- J^sip - 1, q; x, y) 

+{m) - ml)J^siP, q; X, y) - ^J^sip, q;x + 2, y) 

+ ^-f[J^s{P:q-2;x,y) 



-^^s{p,q]x,y + 2) 



Tsip,q;x,y,z,2) 



-2mlTs{p, q-l;x,y)+ mlTs{p, q; x, y) 
^s{p, q-'i-;x, y, z) - Tsip - 1, g; x, y, z) 



+ {ntf - ml)Ts{p, q- x, y, z) - -Ts{P: q]x + 2,y,^ 
-^^s{p, q;x,y + 2, z) - ^J^s{p, q; x,y;z + 2) 



J^sip, g - 2; X, y, z) - 2mlJ^s{p, q-l;x, y, z) 



+mlJ^s{p,q;x,y,z] 
Finally let us notice that we can write, for r > 3 

(kir ^kir' - 4(2 + r'^k')ikir-' sln^ k^ 



(4.13) 



DF{k,mf)P+^ 



DF{k,mf)P+s 

<klY-' . , d 

sm kw 



p + 5-1 "dk^ Dpik, mf)P-^+^ 
Integrating by parts we obtain the recursion relation (to be applied for r > 3) 

^s{p,q;- ■ ■ ,r) 

= 6J^s{p, g; . . . , r - 1) - 8J^s{p, g; . . . , r - 2) - 
4r^^5(p,g-l;...,r-2) + 4r^m^JF5(p,g;...,r-2) + 



(4.14) 
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rw^sip, g - 1; . . . , r - 1) - r^^mlJ^sip, g; . . . , r - 1) + 

, ^ , [-2g-F5(p - 1, g + 1, . . . , r - 2) + ^^^(p - 1, g + 1; . . . , r - 1) + 
p + — 1 I 2 

(2r - 5)J^s{p - 1, g; . . . , r - 3) - ^(r - 2) - 1, g; . . . , r - 2)] (4.15) 

Notice that when this recursion is used for p = 1, terms of the form ^^-^(O, ■ ■ are 
generated: as we aheady noticed in the introduction to the section, this forces us to 
compute J-'s{p, . . .) including terms of order 6 when p < 0. 

Using the previous recursion relations we can now reduce each integral J-'s{p, g; n^, Uy, n^, rit] 
to a sum of the form (setting again mh = rrif = m) 



J^s{p,q]n^,ny,nz,nt) 



p q+k 
r=p—k+l s=q—k 



(4.16) 



where k = + Uy + + rit). It is easy to see, from the structure of the recur- 
sion relations, that, the coefficients ar.s(m, 5) have the properties mentioned at the 
beginning of the section. 

As in the bosonic case we can simplify this expression if we consider the limit 
m 0. Let us first compute the divergent part of the integrals J^s{p,<l) (of course 
we must have p + g > 2). Let us define ADF{k,m) = DF{k,m) — DB{k,m). Then 
let us rewrite 



^s{p,q) 



d^k 
(2^ 



1 



DF{k,m)P+WB{k,m)i 



p+q-2 



DB{k,my+^+' § 



d'^k 



-P 



5\ fADF{k, 




\ I / \ DB{k,m 
ADF{k,my 




{27r)^DB{k,m 



\p+q+l+S 



(4.17) 



It is easy to see that the first term is finite. Thus if we want to compute the diver- 
gences of J^s{p,q) we can limit ourselves to consider the second term which contains 
purely bosonic integrals. Expanding ADF{k,my we see that we need to compute 
the divergent part of purely bosonic integrals Bs{r;nx,ny,nz,nt). If the original in- 
tegral in ( |4.17| ) has p > we need only the ^-independent divergent part, while for 
p < also the terms of order S are needed. The computation of the divergent part 
of Bs{r;nx,ny,nz,nt) can be done in different ways. One possibility is using the re- 
cursion relations of the previous section and the expression for the divergent part 
of Bs{r) reported in ( p. 12 ). One can also attack the problem directly. Indeed if we 



define constants bi(nj.,ny,n;z,nt) by the asymptotic expansion 



n 

i={x,y,z,t} 



Vo(«)) 



-a)' 



E 

i=0 



hi{n^,ny,nz,nt) 



a 



i+2 



(4.18) 
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and q = r — Tlx — Uy — Uz — rit the divergent part of Bs{r; nx,ny,nz,nt) (of course for 
g > 2) is given by 

1 (5 1 \ 

-bq-2inx,ny,nz,nt) ^^^^^^ logm^ M - -logm^ - 5^ - j (4.19) 

From this expression we immediately see that the divergent part of J-'s{p,q) has the 
generic form 

V^^\m){l-5\ogm'^) + 6V^^\m)+ log terms +0(5^) (4.20) 

where T)^^\m) and T)^'^\m) are polynomials in 1/m? whose coefficients are rational 
numbers multiplied by l/vr^ and "log terms" indicates terms which diverge as a power 
of log . 

Exactly as in the bosonic case, the knowledge of the divergent part of J-'s{p, q) can 
be used to simplify ([4.16|) : indeed whenever an integral is multiplied by the infrared 
regulator we can substitute it with its divergent part. Thus we can rewrite ( |4.16|) as 

P g+k 

J^5{p,q;nx,ny,nz,nt) = ^ ^ a,.s(0, 5) ^'^(r, s) + 7^(m, 5) + O(m^) (4.21) 

r=p— fc+l s=q—k 

If p > and g < 0, as only terms with s < can appear in ( |4.16| ), the only ^^(r, s) 
that can contribute to 7l{m, 6) have r > 3. Since ars{m, 6) is finite for S —>■ 0, we see 
that 71(171, 6) is a polynomial in finite for 5 — »• 0. For p > and g > also 

Ts{r,s) with r < can contribute to TZ{m,5). As ars(m, 5) can behave as 1/5 we 
have (see (||20D) 

7^(m,(5) = \ (1 -51ogm2)7^W(m) + n^^\ni) + 0{6) (4.22) 



5 

where 'KP^{'m) and TiS^^m) are polynomials in 1/m^ whose coefficients are rational 
numbers multiplied by l/vr^. 

For p < and g < it is easy to see that TZ{m, 5) = as only finite integrals 
appear in the r.h.s. of ( ^.16] ) while for p < and g > we have 

7^(m,(5) = (1 -51ogm2)7^W(m) + 6n^^\m) + 0{6^) (4.23) 



Thus in all cases the function TZ{m, 5) has the form stated in the introduction to the 
section. 



4.2 Second step: the identity Xi{p,q) 

We will obtain here a new set of recursion relations using the identity Ti{p, g). Ap- 
plying the previous recursion relations we can write each term in ( [4.5| ) as in ( [4.2 1| ) 
obtaining a non trivial relation of the form 

E frsiP, S)^5{r, s) + Usip, g; m, 5) = (4.24) 

r,s 
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where p — 4 < r < p. Following our discussion of the bosonic case we will use ( 4.24 ) 



to obtain new recursion relations. Let us first notice that in ([4.24|) there is only one 
term with r = p — 4. It has s = q + 4 and 

f ( 32(g + l)(g + 2)(g + 3) 

Thus, ii q ^ —1, —2, —3 we can solve ([4 .241) in terms of Tsip — 4, g + 4). Shifting 
p — > p + 4 and g — > g — 4 we can express J-'s{p, g), g 7^ 1, 2, 3 in terms of ^^(r, s) with 
p+l<r<p + 4. We can then use this relation to express all integrals J-'s{p,q), 
p < — 1, g 7^ 1, 2, 3 in terms of J-'s{r, s) with either < r < 3, s arbitrary or 1 < s < 3 
and p < —1. 

Two observations are in order: 

1. A careful analysis of the recursion shows that in the result the coefficients of 
J^s{r, s) with r = 1, 2, 3 are of order 6. This property is very important: indeed 
it guarantees that when substituting these expressions in ( [4.2 1| ) the coefficients 



of J-'s{p, g) with positive p are finite for S ^ 0. This property would not be true 
if we were using the relation to eliminate also J-'s{r, s) with r > 0. 

2. If we are considering Ti{p,q) with g < then the result is expressed only in 
terms of J^i{r^ s) with < r < 3 and s < 0. 

We could also try to use the same identity to obtain recursion relations which express 
Tsip, q) in terms of Tsir, s) with r < p. To do this we should try to solve the identity 
for the J-'s{r, s) with the highest value of r, namely r = p. However in this case there 
are three terms with r = p, namely J-'s^p, g — 2), J-'s{p, g — 1) and J-'s{p, q) and thus 
we cannot obtain a recursion which decreases the value of p. We will thus proceed in 
a different way. We will solve the identity for ^-^(p, g — 2), shifting g g + 2. This 
is always possible as fp,q-2{p,(l]S) = —l- In this way we obtain a recursion relation 
which expresses Ts{p,q) in terms of Ts{r,s) with r < p and J-'s{p,s) with s > q. 
We will use this recursion to eliminate recursively all the integrals with g < —3 and 
p > 4. The choice of stopping at g = — 3 guarantees that only integrals !Fs{r, s) with 
s < are generated. 

A third recursion relation can finally be obtained by solving ( [4.24| ) in terms of 



J-'s{p,q). This is also always possible as fp,q{p,q',S) = 256. In this way we can 
eliminate all integrals J-'s{p, q) with p > 4 and g > 0. 

In conclusion, using the identity Xi(p, g) we can rewrite every Ts{p, g) as 

^&i.P^ = Crs(p, g; 5)J^5(r, s) + 5(p, g; m, 5) (4.26) 

rs 

where in the sum we have either < r < 3, s arbitrary, or r < —1 and s = 1, 2, 3 
or r > 4 and s = —2, —1,0. It is easy to see that the properties mentioned at the 
beginning of the section are satisfied by Crs{p, g; S) and S{p, g; m, 6). 
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4.3 Third step: the identity 12{p,q) 

Let us now obtain a new set of recursion relations which allow to reduce the remaining 
J^s{p,q) in terms of a finite set of integrals. We will use here the second identity 
^2(p, We will discuss separately four different regions: 

1. ? < 0, <p < 3; 

2. g > 0, <p < 3; 

3. g=l,2,3,p<-l; 

4. 5 = -2,-l,0,p>4. 

4.3.1 The strip g<0, 0<p<3 

Let us first consider the integrals J^s{p, l) with q <0 and < p < 3. 

We start from X2(2, q + 4:). We can use the previous relations to obtain an identity 
which involves only jF5(r, s) with < r < 3. If g < — 5 we can solve it for ^5(3,5) 
expressing it in terms of J-'s{r,s) with r < 2 or r = 3 and s > q. Notice that by 
stopping at q — —5 all integrals are expressed in terms of ^^(r, s) with s < 0. 

We want now to obtain a relation for J-'s{2, q) . In this case we start from22(l, Q'+4). 
Then we can use the relations of the first two steps to obtain a recursion relation 
which involves only ^^(r, s), < r < 3. Then we use the previous relation to 
eliminate ^^(3, g — 2), Ts{?>, q — ^) and ^^-5(3, q). At the same time also Ts{2, ? — 2) 
and ^5(2, q — 1) cancel. If g < —5, this relation can then be solved in terms of Ts{2, q), 
the result containing only J\5(r, s) with s < 0. 

In a completely analogous way we can derive recursions for T^O-.q) and ^^^(O, q). 
In the first case we start from X2(0, q + Q). We first apply the step-one and step-two 
relations, then use the previous relations to eliminate ^5(3,5 — 2),J\5(3,g — 1),JF5(3, g) 
and !Fs{2,q). If 5 < — 5 we solve for ^^(l,^). Finally starting from X2{—l,q -\- 6), 
eliminating jF5(3,g — 4), jF5(3,g — 3), .^5(3,5 — 2), J^5(2,g — 2) and ^^(l,^' — 2) we get 
a relation for ^^-^(O, q) vahd for q < —5. 

Using iteratively these four relations we are now able to express every J-^sip,q), 
< p < 3, q < —5, in terms of J-s{r, s) with < r < 3, —4 < s < 0. 

4.3.2 The strip q> 0, <p <3 

Here we want to obtain relations analogous to the previous ones but which decrease 
the value of q. The procedure is identical to the previous one. 

We first consider l2{2,q + 1) and we use the step-one and step-two relations to 
obtain an identity involving only J-s{r, s), < r < 3, then we solve for ^^(3, q). We 
use this relation for q > 0. 

Analogously from X2(l,g + 1) and T2{0,q + 1) we get relations for J^s{2,q) and 
^5(1, q) respectively vahd for 5 > 2 and q > A. Finally we consider X2(— 1, ? + 1): we 
apply the step-one and step-two substitutions, then we use the previous relations to 
eliminate J-s{3, q — 3), J%5(2, q — 2) and J-5{1, q — 1) and finally solve for J-s{0, q). This 
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relation is valid for g > 6. In this way we express all J-s{p, g), < p < 3, g > in 
terms of J-s{r, s), < r < 3, s < 2(3 — r). 

4.3.3 The strip p < -1, g = 1,2,3 

Here we consider X2(p + 3, 1), apply the step-one and step-two relations to obtain an 
identity involving only Ts{r, s) with s < 3, then we solve for J-s{p, 3). This relation 
is valid for p < — 1. In the same way starting from X2(p + 2, 2) we get a relation for 
J^s{p, 2). p < —2. Finally starting from T2{p +1,3), after eliminating Tsip —1,3). wc 
get an identity for J-'s{p, 1), p < —1. Thus all integrals in this region but J-'s{— 1,2) 
can be rewritten in terms of J-'s{r, s), r > 0. 

4.3.4 The strip p> 4, g = -2, -1,0 

Here we start from T2{p — 1,0). We apply the step-one relations and then the step- 
two relations to ehminate Ts{p, —4) and !F5{p, —3), then we solve for ^^(p, —2). This 
relation is valid for all p > 4. Analogously, starting from X2(p — 1, —1), eliminating 
jF^(p, r), —5 < r < —2, we get a relation for Ti,{p. — 1), p > 4. Finally starting from 
X2(p, —2), eliminating J^s{p + l,r), — 6 < r < — 1 and J^s{p,r), — 4 < r < — 1 we get 
a relation for J^s{p, 0), p > 4. In this way all J^s{p, q) in this strip get rewritten in 
terms of Ts{r, s), < r < 3, s < 0. 

4.4 The fourth step: the last relations 

In the preceding step we showed that all integrals can be rewritten in terms of a 
finite number of them. However we have not used the identities Xi(p, q) and X2(p, q) 
for all possible values of p and q. For instance we never used Xi (p, q) for p < 3 and 
q = —1, —2, —3 or p > 4 and q = 0; much larger is the number of cases where the 
second identity has not been used. We thus checked systematically if there were other 
values of p, q for which the two identities were not trivially satisfied, thus providing 
relations which could be used to further decrease the number of independent integrals. 
Using Xi(p, q) with (p, q) getting the values (3, -1), (2, -1), (1, -1), (0, -1), (3, -2), 
(2, -2), (1, -2), (1, -3) and X2(2, 0) , X2(3, -2) we obtain relations for J^5{3, [0, -1]), 
J^s{2, [-3, -4]), J^s{^, [-3, -4]) and J^5(0, [-1, -2, -3, -4]). These relations arc re- 
ported in the appendix. Notice that each integral is expressed in terms of J-'s{r, s) 
with s < 0. We stress that the identities we used to derive the relations were chosen 
arbitrarily; other choices could have been equally used. However, once these rela- 
tions have been computed we have found that Xi(p, g), (resp. X2(p, g)), is identically 
satisfied for all values of g and for all p < (resp. p < 0). 

In a completely analogous way we have found that the identities Xi(4,0) and 
X2(p,g) with (p,g) getting the values (3,2), (3,3), (1,2), (1,3), (0,4), (0,5) are 
not yet satisfied. We have used them to get relations for Ts{2, [1,2]), ^^(1, [3,4]), 
J-'s{0, [4,5,6]). They are reported in the appendix. 

Collecting everything together wc get 

^s{P: q; n^^: %, nnt) = ^ drs{5)Ts{r, s) + T{m, 6) (4.27) 
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J^(1,0) 


0.08539036359532067914 


^(1,- 


-1) 


0.46936331002699614475 


^(1,-2) 


3.39456907367713000586 




-1) 


0.05188019503901136636 


^(2, -2) 


0.23874773756341478520 


m- 


-2) 


0.03447644143803223145 


^(3, -3) 


0.13202727122781293085 




-4) 


0.75167199030295682254 




- 0.01849765846791657356 






0.00376636333661866811 


Y2 


0.00265395729487879354 


Y; 




0.00022751540615147107 



Table 2: Numerical values of the constants appearing in the fermionic integrals. 



where drs{5) and T(m, 5) have the form explained at the beginning of the section. 

Now, if p > and g < the second sum extends over ten values of (r, s), i.e. 
(0, 0), (1, t), (2, t) and (3, t + 2) with < t < 2. Let us consider the limit 5-^0. The 
polynomial T(p, g; m, 5) is finite in this limit and the same is true for all drs{p, q] S) 
with r > 0. Thus the only coefficient which could have a. 1/6 divergence is doo{p, q; S). 
However the result is finite for 6^0 and thus also this coefficient is finite in this 
limit. We can thus set simply 6 = and use the fact that ^^-^(O, 0) = 1 + 0{6) to get 

J^{p, q; n^, Uy, nut) = ^ drs{p, q; 0) J'(r, s) + doo{p, q; 0) + T{p, q; m, 0) (4.28) 

r>0;s 

Finally let us consider the limit m — ^ 0. In this limit all J-'{r, s) appearing in the 
previous sum are finite except JF(2, 0) which we write as 

^(2,0)= -^(logm2 + 7E-Fo) +Fo (4.29) 

where Iq is a numerical constant. Notice that if J^{p,q;nx,ny,nnt) is finite, since 
JF(2, 0) is the only term which can contain logm^, d2fi{p,q]0) = 0, i.e. the result 
depends only on the eight finite integrals. Numerical values are reported in Table 

Let us now consider the case p > and g > 0. In this case in the sum on the r.h.s 
of (g^) we can also have J^s{-h 2), ^"5(0, 1), ^^^(0, 2), ^^^(0, 3), ^"^(1, 1) and ^^^(1, 2); 
T(p, q; m, 6) has the form 

\t^^\p, q; m)(l - (51ogm2) + T^^\p, q; m) + 0(5) (4.30) 


where T^^\p,q;m) and T^'^\p,q;m) are polynomials in 1/m^. Let us consider the 
limit 5 — > 0. The only coefficients that may behave as 1/5 are those with r < 0. 
Writing in this case 

drsiS) = Idi"^ + d^^J + 0{6) (4.31) 


the cancellation of the 1/6 terms gives the equation 

4)^(3) + rfSi3(2) + (rfff + 2d%)Bil) + + T(i)(m) = (4.32) 
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where we have used ^^-^(O, q) = B{q) + 0{6) and 

^•^(-l, 2) = B{1) - B{2- 2) + 1/4 + 0{5) = 2B{1) + 0{5) 
In all cases we have found that ( |4.32| ) is satisfied in a trivial way, i.e. 

"03 ~ "02 — "OO ~ U 

4^ + '^d% = 



(4.33) 







(4.34) 



This is not completely surprising: indeed cancellation of the terms divergent for 



m 



requires d^ 



(1) 

03 



— 8(io2'* and 



m] 



(1) 

03 



327r^m^ 



+ 



(4.35) 



where t is a rational number. Then (|4.32| ) becomes 



(1) 

03 



Z, + {dS + 2d'LU)Zo + d\ 



256 



(1) 



00 



d'i 



(1) 

03 



1024 



+ 



13 
1536' 



"03 







(4.36) 



Thus, if for some p,q,nx,ny,nz,nt, the relations ( [4.34|) were not satisfied, one would 
get an equation which would allow to eliminate either Zo or Zi. 

Because of ( |4.34|) , in the limit 5 — > 0, we can then substitute in (f4.27|) ^^-^(O, 3) = 
B{3), J^5(0,2) = B{2), J^s{0,0) = 1 and rewrite 

doi^5(0, 1) +ci_i,2^5(-l, 2) = idoi (2^5(0, 1) -^5(-l, 2)) + {2d.i,2 + doi)B{l)+Oi5) 

(4.37) 

showing that ^^-^(0, 1) and J-'s{—l, 2) appear in the result only in the fixed combination 
2J^s{0,l) — J^s{—^,2). In conclusion, for g > the result beside the nine integrals 
which appear for g < contains also the bosonic constants Zq, Zi, Fq — •je and the 
integrals ^-"(1, 1), ^(1, 2) and lims^o{2 J^s{0, 1) — ^j(— 1, 2))/S. As in the bosonic case, 
instead of these three quantities we have parameterized our results in terms of three 
infrared-finite integrals. We introduce: 



= -^(1,1; 1,1,1) 

Y2 = -^^(1,1; 1,1, 1,1) 
16 

>3 = -^^(1,2; 1,1,1) 
16 



(4.38) 



Their numerical values are reported in Table ^ The relation with the original integrals 
is: 



lim^(2^5(0,l)-^5(-l,2)) = ^-12^1-3^0 + 2^(1,0) 

5— >0 d 4 



(4.39) 



+— -2) + — J^fl, -1) + — 0) + — J^(2, -2) - — J^(2, -1) 
768 ^' ^ 192 ^ ' ^ 192 ^ ' ^ 768 ^ ' ^ 48 (4.40) 

^/ ^ 1 19 307 3 19 19 1 , 
TT/^i 2^ = I u V V -I V -I V 

^ ' ^ 327r2m2 12288 184327r2 64 ° 256 ^ 1024 ^ 8 ^ 
1 19 53 187 

1 73 

-^J^(3,-2) (4.41) 

2304 ^ ' ^ ^ ^ 

As in the bosonic case, let us notice that logm^ appears always in the fixed combi- 
nation (logm^ + 'Je — Fq). Thus in finite integrals Fq — does not appear. 

To conclude this section we want to add a few remarks on the numerical evaluation 
of the constants appearing in Table 0. A direct evaluation of the integrals does not 
provide accurate results: we have thus used a different procedure inspired by the work 
of . To evaluate the constants q) with g < we have considered the integrals 
Jq = JF(1, —q) with 6 < g < 13. To compute them we have first calculated the sums 

where k runs over the points k = {2n/L){ni + ^,n2 + |,n3 + |, ^4 + |), < rij < L 
for various values of L between 50 and 100. Then we tried to extrapolate Jg^L using 
the form 



4^ = (,1 + j (4.43) 

In all cases the correction turned out to be completely negligible. In practice Jq could 
be determined with a relative error < 10^^^. Then each Jq was expressed in terms of 
the basic constants. We obtained in this way 8 equations which were solved giving 
the results of Table 0. 

Analogously to compute Yq, Yi, Y2 and Ys, we computed numerically JF(1,1,8), 
J-'{2, 1, 9), J-'{3, 1, 10) and ^^-"(5, 2, 11) and then solved the corresponding equations. 



5 Integrals in other infrared regularizations 

In the preceding two sections we have discussed the computation of bosonic and mixed 
bosonic-fermionic integrals using as infrared regulator a mass m. Here we want to 
discuss other types of infrared regularization: first we will consider the dimensional 
regularization ||^ and then we will consider mixed bosonic-fermionic integrals with 
the exact Wilson-fermion propagator (|4.3| ). 
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5.1 Dimensional regularization 

In this case we consider the integrals M 



T^^{p,q]n-,,ny,n^,nt) 



(5.1) 
(5.2) 



It is easy to see that the basic recursion relations (p.4|) , (pl7|), ([4.10|) , ( |4.13|) , ( [4.15|) 
can be easily generalized to dimensionally-regularized integrals. The relation ( p.l6|) 
and ( |4.5|) , ([4.6|) are instead intrinsically four-dimensional identities and for this reason 
we have made the computation using a mass as a regulator. Now, we will show how 



to compute ( |5.1| ) and ( |5.2|) from their mass-regularized counterparts. 

Let us begin with the bosonic case, considering B{p] rix, Uy, n^, rit). If q = p — Ux — 
Tiy — rix — rit < 2 the integral is finite and thus independent of the infrared regulator. 
For g > 2 let us rewrite 



B{p;ncc,ny,n^,nt) 

Then define 

F{p;nx,ny,n^,nt) 



1 



2PT{p) Jo 



2PT{p) 
1 



daa^-'e-''^"/^ U 

i={x,y,z,t} 



daaP-^ n 

i={x,y,z,t} - 



(e-"/o(a)) 



(5.3) 



d" 



da^- 



2PT{p) Ji 



da aP-^ \ n 

J={x,y,z,t} 



(e-%(a)) 



d"" 



(-2)"-— ^(e-%(a)) 



a 



E 

i=0 



a 



i+2 



(5.4) 



where hiin^, ny,nz,nt) are defined by ([4.18|) . Then for — we have 



B{p; n^, Uy, n^, rit) = F{p] n^, Uy, n^, rit) 



1 



9-1 



1 



^h~2{nx,ny,n:,,nt) 



T{q-i) 



1 1 



^2g-2i 21 q-i 



Tfl 

^g-^^^^bg-2{nx,ny,n„nt) ( log — + 7ij 



(5.5) 



Analogously we have 



S^^(p; nx,ny,n,,nt) 



1 /""^ 



2PT{p) 



J daaP-^ [e-'^Ioia 



(i-4 



n 

i={x,y,z,t} 



-21 



da"' 



[e-%(a)) 
(5.6) 
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so that for e = 4 — c? ^ we have 
2 1 



bq-2{n:c, Uy, U^, Ut) (1 + eCq-2{n:c, U^, Ut)) + F{p; Tlx, Uy, n^, Tit) (5.7) 



e 23r(p) 

where Ci{nx, riy, n^, rit) are defined by 



n 

i={x,y,z,t} 



log(e °'Io(a)^/a) 



«=0 



Q! 



(5.8) 



Then, by comparison, we have 

-Y^Z_,h-2[n:c,ny,n^,nt) y 



29r(p) 



2i ^2g-2i 29 (g - i)y 

1 ( rn? 2 \ 

bg-2{nx, Uy, n^, nt) - log — - 7b H h 2cg_2(n^, n^, rit) (5.9) 



A simphfication occurs if the integral is logarithmically infrared divergent, i.e. if 
q = 2. In this case, as co{nx, riy, n^, rit) = | log27r, we have 



B^^{p;nx,ny,n^,nt) = B{p;nx,ny,n^,nt) 



4r(p) 



bo{nx,ny,nz,nt) 



log m — 7e H h log 47r 

e 



(5.10) 



Thus, for these integrals, we can use a very simple recipe to go from the mass regu- 
larization to the dimensional one: simply substitute in each integral log + 7^ with 

2/e + log 47r. 

Let us now consider the fermionic case. Again we should consider only the case 
Q = p + q — Ux — Uy — Uz — Ut > 2. Then let us rewrite 



•^^^(P, nx, Uy, nt)= j 



^2n4.2»,fc2n4.2n. Q_^2 / N / AD^{k^ 



i.2nx i.2ny i.2nz i.'int 
ftj, Avj. Avj 



Dsik, 0)^'+9 1^ 



DB{k, 0) 



Q-2 



+ E 

1=0 



I J 



f ^'^^ £.2n, L2ny pn, i2nt ^Dpik^^Y , . 



where ADp{k,m) = Dp{k,m) — DB{k,m). The first integral in the is clearly finite. 
We can thus set e = 0. Then it is easy to see that we can add everywhere a mass, i.e. 
substitute £'b(/c,0) with DB{k,m) and analogously for £)i?(/c,0), without changing 
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its value in the limit m — 0. Thus we get finally 



J^^^ip, q; n^, Uy, n^, rit) = J^ip, q; n^, Uy, n^, rit) 



+ V I ^ 1 i /" pn,l2nyl2n,f 2nt 

(2^)' ' Ds(A;,0)f+'?+' 
- /7^^r^^r^^r^^^^^^%4| (5.12) 

In the last term purely bosonic integrals are involved and we have already discussed 
how to compute the difference between their value in the two regularizations. Notice 
that, as in the bosonic case, logarithmically divergent integrals can be dealt with 
easily: simply substitute (2/e + log47r) to (logm^ + 'Je)- 

5.2 Massive Wilson-fermion-propagator integrals 

Here we want to consider integrals of the form 

T{p,q;n^,ny,n^,nt} = — , ^ 5.13 



where Dp{k, m) is defined in ( ^4.3| ). When Q = p + q — — Uy — — Ut > 2 the 
integrals diverge for m — 0. We will now relate them to JF. Indeed we can rewrite 



J^{p, q; rir,, Uy, n^, rit) = ^{p, q; n^, riy, n^, rit) 



m 



where ADp{k, m) = Dp(k, m) — Dp(k, m) = mk"^. 

Let us notice that if the integral is logarithmically divergent {Q = 2), for m ^ 
we have T = T . For integrals with Q > 2 we see from the explicit expression ( ^.14|) 



the the divergent part is now a polynomial in 1/m instead of For this reason 

the expression for T are in general more cumbersome than those involving T and this 



is why we have studied integrals with (|4.2|) instead of 



6 Applications 

In this section, as an application of our method we will give a few examples. In the first 
subsection we will report analytic expressions for various renormalization constants 
whose value is reported in the literature only in numerical form or is expressed in 
terms of cumbersome integrals. Then we will briefiy discuss the computation of the 
fermionic propagator. 
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6.1 Analytic expressions 
6.1.1 Fermionic self-energy 

We want to give here the expression for the fermionic self-energy. The first computa- 
tion for the Wilson action in Feynman gauge was given in (]Tl[ and it was subsequently 
corrected in [|1^ ]. The fermionic self-energy at one loop has the generic form 

E^^^^(p2,m2) = g^^f--^(5m + i^T.^{p\m^) + mT.2{p\m^)) (6.1) 
For Tw = 1, in Feynman gauge, in terms of our basic integrals we have 



5m = -Zo - 2J^(1,0) ^ -0.3257141174 (6.2) 

l/„^/00 OOn 1^ 1 1 ,^ 1^ 

1 



Ei(plm2) = -(2G(p'^a^,m^a^) +-yE - Fo) + -Zo + — ^—-Yo + -Yi 

' ^ 167r2^ ' ^ ' "^ 8 192 327r2 "^4 



\^(2,-2) + ^.F(2,-l) 



768 ' ' ' 48 
~ -^G(pV, m^a^) + 0.0877213749 (6.3) 

E2(p^m2) = -L(F(pV,mV)+7i.-Fo) + ^-^-4Fo+ ^-1-^1^2 

^J^(l, -2) - —m, -1) - — J^(l, 0) - — J^(2, -2) 

192 ^ ' ^ 48 ^ ' ^ 48 ^ ' ^ 192 ^ ' ^ 

49 1 
+ —T{2, -1) ^ —F{p^a'^, m^c?) + 0.0120318529 (6.4) 
12 47?^ 

where 

F (p'^ , a^) = / (ix log[(l — x)(p^x + m^)a^] (6-5) 

JO 

Gip^a^ ,rn?o?) = / a; log[(l — x)(p^x + m^)a^] (6.6) 
6.1.2 Gluonic self-energy 

Let us now consider the gluonic self-energy which was firstly computed in ||. 
The contribution of the fermions, for ry/ = 1, can be expressed in terms of our basic 
integrals as 



Nf 



127r 



^ '6 if (p^a^m^a^) +7i5 -log 47r) +L 



(6.7) 



where Nf is the number of fermions which are in the fundamental representation of 
SU{N), m the fermion mass (for simplicity we assume all fermions to have the same 

^ Notice however that formula (3.15) in jlj] contains a misprint: the correct result is given in 
formula (10b) of O. 
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mass) , 



-1) - j^(l,0) + A^(2, -2) + H^(2, -i; 



and 



H{p^o? ,m^a^) = I dx x{l — x) \og\x{l — x)p^ + 
Jo 



2 21 

m a 



(6.8) 
(6.9) 



Numerically L ^ 0.0031048512 in agreement with 
For p — »• we have 



n^(p, 



m] 



^/.2. 2 



9 {P ^t^u - Pt^Pu 



127r2 



logm^ - 0.013391999 + 0{p^ 



while for m = we have 



n^(p,o)= -^/(A 



Pf^Pu) 



logp2 - 0.027464385 



(6.10) 



(6.11) 



6.1.3 Renormalization constants for bilinear fermion operators with the 
clover action 

We want to give here the renormalization constants of bilinear fermion operators with 
the clover action 0, |1^ , using the explicit expressions in terms of lattice integrals of 
| T5| . We define local operators 

^LATT^ioc^^^ = ^{x)Ti;{x) (6.12) 

and improved operators 
^LATT.m,^^-^ = ^(x)r^(x) + ^ E (D^ii{x)^^Ti;{x) - ii{x)T^,D,i;{x)) (6.13) 



D^il){x) = - \Ui,{x)iIj{x + /i) - U^{x - ij,)iIj{x - n) 



(6.14) 



where F is a Dirac matrix and 

1 
2 

For each operator we compute a finite renormalization constant Z such that 

(/|0^«^^(x)K) = Z{f\0^^^^{x)\^ (6.15) 



where / and i are arbitrary external states. In the continuum we adopt the MS- 
scheme with scale = 1/a. We write at one loop 



1 + ^? 



(6.16) 



Expressions for AZ, in Feynman gauge, are reported in |TB| in terms of quite compli- 
cated integrals. The expression for ry/ = 1 in terms of our basic integrals is reported 
in Table |^ and Table The final numerical values are in agreement with those of 
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A ^75 7m 




1 




31 


7 


3 


11 


23 




48 


48 


32 


32 


144 


l/7r2 




223 


139 


19 


5 


73 




128 


128 


64 


8 


384 




7ij)/47r2 


-3 


-3 








1 


Yo 




-12 


-12 


-2 


-2 


4 
3 


Yi 




-18 


-6 


23 
2 


35 
2 


46 
3 


Y2 




1 

4 


1 
4 


3 

8 


3 

8 


7 
12 


Y3 




-48 


-48 


-48 


-48 


-48 


Zo 




6 


2 


3 


5 


4 


Zi 




3 


3 


3 


3 


3 




4 


4 


4 


4 


4 




-2) 


5 

96 


7 

96 


15 
128 


7 

128 


1 
9 


-^(1,- 


-1) 


391 
192 


139 
192 


39 
32 


15 

8 


923 
576 


^(1,0) 


53 
8 


99 
8 


1031 
96 


755 
96 


665 
72 


m- 


-2) 


65 
192 


271 
192 


235 
128 


123 

128 


971 
576 




-1) 


407 
8 


117 

8 


631 
24 


155 
24 


1387 
72 




-3) 


77 
32 


105 
32 


7 
2 


49 
16 


329 
96 




-2) 


1057 
16 


203 
16 


259 

8 


7 


413 

16 


Total 




-0.4891266 


-0.5669576 


-0.3882898 


-0.3493743 


-0.2819883 


Ref. 1 


15 


1 


-0.495 


-0.573 


-0.388 


-0.349 


-0.280 



Table 3: Results for AZy for the local operators (|6.12|) : we report here the coefficients 
of the various constant appearing in the result. "Total" is the numerical value of AZp. 
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AZid 


75 




151 IJ. 


C> jj,i/ 


1 




151 


7 


19 


41 


11 




48 


48 


24 


48 


24 


l/7r2 

/ 




279 

128 


139 
128 


3 
16 


47 
64 


101 
384 


(Fo- 


7£)/47r2 


-3 


-3 








1 


Yo 




-12 


-12 


-2 


-2 


4 
3 


Yi 




-18 


-6 


1 
2 


_ n 
2 


2 
3 


Y2 




_i 

4 


_ 1 

4 


3 

8 


3 

8 


J_ 

12 


Y3 




-48 


-48 


-48 


-48 


-48 


Zo 




6 


2 


-1 


1 


_4 

3 


Zi 




_3 


_3 


_3 


_3 


_3 




4 


4 


4 


4 


4 




-2) 


4 

Q 




_ n 


77 


23 




-1 


) 


329 


139 


77 


113 


11 




192 


192 


384 


384 


32 


-^(1,0) 


157 

8 


_99 
8 


1373 
96 


163 
96 


173 

18 


m- 


-2 


) 


4321 


271 


2737 


1855 


323 




192 


192 


384 


384 


64 




-1) 


1823 
8 


117 

8 


1693 
24 


1217 
24 


3511 
72 




-4) 


7 
32 





7 

128 


7 
128 


7 
192 


m- 


-3) 


1585 
32 


105 
32 


1055 
64 


635 
64 


145 
12 


m- 


-2) 


4713 
16 


203 
16 


179 

2 


513 

8 


3067 
48 


Total 




-0.2634473 


-0.5669576 


-0.2493438 


-0.0975887 


-0.0591138 


Ref. 1 


15 


1 


-0.269 


-0.573 


-0.249 


-0.0973 


-0.0570 



Table 4: Results for AZ^ for the improved operators ( |6.13| ): we report here the 
coefficients of the various constant appearing in the result. "Total" is the numerical 
value of AZp. 
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6.1.4 Renormalization constants for fermionic energy-momentum tensor 

We want now to compute the renormalization constants for the dimension-four op- 
erators which appear in the first moment of the deep-inelastic-scattering structure 
functions and in the definition of the energy-momentum tensor. We will consider two 
operators whose continuum formal limit is 



^ ip-ff,D^ip - D^ip^^ip + in^u) (6.17) 



For their explicit definition in terms of lattice operators we refer to |JT6[, formula 



(5.22) and [||], formulae (3.14)/(3.16). We want now to compute new operators so 
that 

(/lag'^^^^K) = (/loS'^^^^K) (6.19) 



for arbitrary states / and i. In the continuum we adopt the MS'-scheme with scale 
jjL = 1/a. To define O^^l'^^"^^ we must consider all possible mixings with operators 
of dimension less than or equal to four. Here we will restrict our attention to the 
gluonic sector. We will write at one loop 



^^LATT ^ 0W + (g_20) 



where 



1 X LATT /I X LATT 



LATT 



The superscript LATT indicates that we use here some lattice operator with the 
given continuum limit: the explicit discretization is however irrelevant at one loop. 
For rw = 1 the constants are reported in table ||. For the gluonic operator ( |6.18| ) 
we write A^'^ = Af^^ + NfA\^'^\ The constants A\^^^ have been computed in see 
formula (6.11)0. As for the constants Af^\ we have A^^-^^ = 0, A^^-^^ = Af^\ The 
explicit expression of Af for rw = 1 is given in table 

We can also easily compute the renormalization constants for the operators which 
are the trace of O^^J and O^^J, i.e. for 

^^'^ = ^E(V^7m^m^-^m^7m^) (6.22) 
= F:,F^, (6.23) 



^The calculation in considers minimal subtraction in the continuum. If one considers the 
MS'-scheme, one should use the formulae of sect. 6, with Y = {Fq — ^/e)/'^^- 
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A« 


a(i) 

^3 


aS'^) 


1 


2 


1 


11 


1 




27 


108 


108 


9 


l/7r2 


11 
144 


1 
288 


19 
96 





(i^o-7E)/47r' 


1 

3 








1 
3 




4 
3 








4 
3 


^(1,-2) 


1 

54 


1 

432 


11 
432 


1 

36 


^(1,-1) 


2 

27 


25 
432 


47 
432 


1 

18 


^(1,0) 


13 
9 


7 

36 


137 
36 


7 
6 


^(2, -2) 


1 

54 


97 
216 


43 
216 


5 

24 


^(2,-1) 


62 
9 


13 
6 


359 
18 


2 
3 


^(3, -4) 





13 
288 


13 

288 





J^(3,-3) 





1 

8 


17 
72 





^(3, -2) 


149 

18 


113 

36 


791 

36 





Total 


0.00826199 


-0.01058036 


-0.00963232 


0.01339200 



Table 5: Results for A^'-*: we report here the coefficients of the various constant ap- 
pearing in the result. "Total" is the numerical value. All constants must be multiplied 
by Tf defined by Tr T^T^ — Tf5'^^. For fermions in the fundamental representation 
of SU{N) we have Tf^\. 
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Considering again only the gluonic sectorQ we can write 

0{3)_^^2^(3)(^(4) 



0(3) 

aw 



We haveQ 



A(3) 
A(4) 



247r2 



lliV 
967r2 



(6.24) 
(6.25) 



(6.26) 
(6.27) 



I for fermions in the fundamental 



where Tf is defined by Tr T"T^ 

representation of SU{N)). The last numbers in A^^) and A^^) are related to the 
anomaly of the energy-momentum tensor. Indeed define 



(6.28) 



where the first sum is over the Nf fermion species. Then 



T, 



MM 



niN 

\96~ 



247r2 



NfTf ) F 



(6.29) 



in agreement with p^ . 

Finally we want to compare our results with those of Capitani and Rossi |17]. An 
easy calculation gives: 





-''^ N 


4 
3 


BL = 




--Tf 

9 ^ 


Bqg = 




9 ^ 



2tt^ 

-17.778285 + — 
« -2.16850094 (2Tj) 
0.20789614 (2T/) 



(6.30) 
(6.31) 
(6.32) 



As already noticed in [T^ our final result for Bgg agrees with theirs^. We are also 
in perfect agreement for B^g, while we differ for Bgg which is reported in [|T^ to be 



B, 
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0.019 for T 



f 



^ As before there are additional mixings with dimension-three and dimension-four fermionic 
operators. For a numerical evaluation of these mixings see [l^ . 

^The additional terms which appear in A^^^ and A^^^ are due to the fact that in dimensional 
regularization (5^.iV(0^,) ^ 7V(<5^,0^,) 1^. 

^There are however some misprints in Tables 10 and 11 of in Table 10 the contribution 
of "Sails" is -85/(1447r2) - 7/(247r2)(2/e + log47r - 7^) and "Total / - J" is {-1 / {Un^)) , while 
in Table 11 in "Sails" -7/(9tt'^) should be replaced by -S/{16t:'^) and in "Total J", -ll/(187r2) 
should be — l/(487r^). We are also in disagreement with their Table 2, where we get Bgg = —4/3, 
B^g = —20Tf/9, Bgg = —8Tf/9. However a recomputation confirms our results [ p3| . 

^'^ This discrepancy has been recently understood and the new result is in agreement with ours [p3[ . 
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We want also to correct here the results which appeared in [24|. Indeed the 
constants Z5, Zq and Z-j appearing in Table 1 should be Z5 = 1 + (7^1. 02165, Z^ = 
—^f^O. 65205 and Z7 = (7^0. 25034. The different result for Z5 and Zj was due to a 
numerical mistake in the evaluation of the photon contribution (the exact result is in 
pi], sec. 5). There are also minor misprints: the sign on in (16), (17) and (18) of 



24 1 should be "minus". An analogous sign should be changed in section 7 of ||16|| . 



6.2 Calculation of the propagators 

Recently a very efficient numerical method to evaluate higher-loop integrals has been 
presented in [Q. An essential ingredient in the method is the exact calculation of the 
free bosonic propagator. In |^ an algorithm was introduced which allowed to express 
its values in x-space in terms of the values at x = (0, 0, 0, 0) and x = (1, 1, 0, 0). Here 
we will discuss how to use our method to obtain an algorithm for the bosonic and the 
Wilson-fermion propagator. 
In general we will consider 

/o Yi n fi \^ (6-33) 

(27r)* DByk, m)P 

77^-- (6.34) 

The computation of these two quantities is in principle straightforward: indeed, us- 
ing the symmetry k^ — /c^ we can rewrite e^^^ as O/x cos(A;^x^) and then express 
cos(/c polynomial in k"^. In this way Gb{p', x) and Gf{p, x) are expressed^ 

in terms of B and and thus, for m — > 0, using the results of the previous Sections, 
we can express Gb{p',x) in terms of Zq and Zi (and Fq if p > 2) and Gp{p,q;x) in 
terms of the eight constants JF(1, [-2, -1, 0]), T{2, [-1, -2]) and JF(3, [-4, -3, -2]) 
if g < and p = 1, to which we must add J^{2, 0) if g < 0, p > 2, and Iq; ^i; ^2, ^3, 
Zq and Zi'ii q> 0. Of course, in a completely equivalent way, in the bosonic case we 
can use instead of Zq and Zi the values of 6*^(1; x) at two different values of x. For 
instance, as in @], we could choose 

Gb(1;(0,0,0,0)) = Zo (6.35) 
Gb(1;(1,1,0,0)) = Zo + Zi--^ (6.36) 

Analogously in the fermionic case we could express Gf{p, Q', x) for g < in terms of the 
values of Gf{1,0;x) at eight different points: a possible choice, with all the points 
in a hypercube of side two, is given by (0,0,0,0), (1,0,0,0), (2,0,0,0), (2,1,0,0), 
(2,2,0,0), (2,2,1,0), (2,2,2,0) an d (1, 1, 1, 1). 

^^Notice that also the reverse is true: the integrals B and T can be rewritten in terms of and Gp 
using the identity (A:^)™ = (— l)™(e*'^/^ — Q-^k/2^2m^ Loosely speaking Gb and B (and analogously 
Gf and J-) are two different "basis" in the space of bosonic (respectively mixed bosonic-fermionic) 
lattice integrals. 
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From a practical point of view, one can simplify the algorithm by implementing the 
"integration-by-part" recursion relations (|3.7|) and ( [4.15| ) directly on the propagators. 
To apply our technique we will generalize ( |6.33| ) and ( |6.34| ) by considering 

(6-38) 

For these two quantities it is very easy to obtain recursion relations. We obtain in 
the two cases: 

gB{p;x + fi) = gB{p;x-fi) — -gB{p-l;x) (6.39) 

p + - I 

gF{p,q;x + 2fj,) = gF{p,q;x - 2fi) - r^^ {gF{p,q ~ l;x + fi) - gF{p,q - l;x - fi)) 

+ {m^r^ - 2m) (^f(p, q;x + fi) - gF^p, q]X- fi)) 
2 

- p_^^^ [Xf^SF{p-l,q]x) 

+q {Gf{p - 1, g + 1; X + /i) - gF{p - 1, g + 1; X - /i))] (6.40) 

where is a lattice unit vector. Using these two relations it is possible to express any 
element ^b(p; x) in terms of gsip'] x') where x' is an element of the unit hypercube 
(i.e. x'fj, = or 1). Analogously in the fermionic case we express every gF{p,q',x) in 
terms of gF{p',q';x') with x' belonging to the hypercube of side two, i.e. x'^ = 0, 1 
and 2. These last quantities can then be easily expressed in terms of Bs or JF^ and 
then the procedure we have been presented in the previous Sections can be applied. 
In the bosonic case, the algorithm we have described is less efficient than the one 
introduced by [^: indeed for p = 1, starting from ( |6.39|) , Vohwinkel obtains a 
recursion relation which involves only terms with p = 1 and which thus avoids the 
necessity of introducing 6 and nonpositive values of p. In the fermionic case however 
we have not been able to implement the same trick. 

Finally let us notice that once ( |6.34| ) has been computed one can easily obtain the 
fermion propagator 

(27r)^ DF{k,m) 

as 

^f{x) = -W^^iGF{l,0;x + fi)-GFil,0;x-fi))- 

^ H ((^^(1, 0; X + /i) + Gf(1, 0; X - ^u)) + [^rw + m)Gir(l, 0; x) 

(6.42) 

To conclude the section let us discuss the numerical evaluation of the fermionic 
propagator. As in the bosonic case 0], the expressions for Gf{p, q] x) in terms of our 
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basic constants are numerically unstable for \x\ — >■ oo: a numerical error in the basic 
constants gets amplified in Gp{p, q; x) as — >• oo. This problem has a standard way 
out: if the expressions are unstable going outward from the origin, they will be stable 
in the opposite direction: thus, if we want to compute the propagator for |x| < d, 
for some fixed d, we choose eight points with |x| ~ d (say yi,...,ys) and then we 
express the propagator for |x| < d in terms of GF(l,0;?/j), i = 1,...,8. The new 
expressions are numerically stable: the numerical error on G'ir(l, 0; gets reduced 
when we compute the propagator for |x| — > 0. As noticed in ^ the instability of 
the recursion can also be used to provide precise estimates for the basic constants. 
We have thus used this method to obtain an independent numerical estimate of the 
first eight constants of table |^ considering the set of eight points X*^") = {{n, [0 — 
3],0,0),(n + 1, [0 — 3], 0,0)}. If one estimates the propagator G'i7(l,0;x) at X^"^ 
with an accuracy of 1%, one reproduces the results of Table § with an accuracy of 
approximately 10^^ (respectively 10^^) with n = 7 {n = 9 resp.). 

It is also possible to apply the same procedure to 6*^(1, g; a;) with g < 0. The 
main advantage is that, using larger negative values of q, one can obtain more precise 
estimates of (7^(1, g; x) at the set of points X^"'\ With q = —2 and n = 7, using 
the method presented at the end of section (i.e. evaluating the integrals by 
computing discrete sums of the form ( |4.42| ) with L = 50 - 100), we obtain the 
values of G'ir(l,— 2;x) at the points X*^"-* with a relative precision of 10~^ and the 
final estimates of the basic constants with a precision of 10~^. If one needs better 
estimates one can obtain them by simply increasing both n and —q. The precision 
can in this way be increased at will. Using n = 26 and q = —3 we checked the first 
eight constants of table 0. 



Acknowledgments 

We want to thank Pietro Menotti for collaborating in the early stages of this work. 
We also thank Stefano Capitani and Gian Carlo Rossi for useful correspondence. 



A Final relations 



In this appendix we report the relations which have been found in the fourth step of 
our reduction procedure for the fermionic integrals. The relations for J^s{p, q), q ^ 
are 



jr^(0,-4) = 7336 ^"5(0,0) + 5 
. 2538989 ^, 



910433 



1447r2 
501607 



1024087 325123 _ ,^ ^, 
H 7-r. — -^(5(1, —2) 



36 



144 



72 ^ ' ' 6 

1919221 , , 7343317 
+ 77— -4) + ^^(3, -3) - 

9029 12392 13138 



3075865 



48 

^•^(0,-3) = 704J?^5(0,0) + 5 



97r2 



27 



27 



36 

20889553 



9 



-^5(3,-2) 



-^5(2,-1) 
(A.l) 
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^ll^Ml, -1) - ^^.(1. 0) - ^^,(2. -2) + ?H|i5^,(2. -1) 



369580^,^ 3903064^,^ ^/ 
+ 1190 ^5(3, -4) + ^^(3, -3) ^5(3, -2) 



(A.2) 



J^siO, -2)=72jr,(0, 0) + S 



3« «3 T^^,,_,) 



47r2 3 12 ' ' 6 
8173 

- 678 j^5(l, 0) - -^-^5(2, -2) + 8340 j^5(2, -1) + ^^^(3, -4) 
+5265 :r5(3, -3) - 15322 J^s{^, -2)] (A.3) 
15 



^5(0,-1) =8^5(0,0) + 5 



- 4^-5(1, -2) + 25 J^5(l, -1) - 228 J^5(l, 0) 



27r2 

-268 J^5(2, -2) + 1608 J^5(2, -1) + 630 ^^^(3, -3) - 2436JP'5(3, -2)] (A.4) 
.^ 1606 6088 2962^,^ 3548^,, 

, , 100312 , , 282272 , , 1366 , 
+ 14992^5(1, 0) + —^^,(2, -2) 3— -^^(2, -1) + ^^6(3, -4) 

86456 _ , , 388304 _ , , , , , 
^^5(3, -3) + ^^(3, -2) (A.5) 

Ul, -3) = ^ + ^ - f ^5(1, -2) + -1) - 528 0) 

X 756 

—J^5{2, -2) + 3696 J^5(2, -1) + 1366 J^s{3, -3) - 5368 J^5(3, -2) (A.6) 

O 

^,(2, -4) = ^ - ^ + :^^,(1, -2) - ^^,(1, -1) - 264^,(1, 0) 

52 

+— j^5(2, -2) + 1200^5(2, -1) - ^5(3, -4) +4^5(3, -3) - 1304^5(3, -2) (A.7) 

^5(2, -3) = +5^5(1, -1) + 36^5(1, 0) - 4^5(2, -2) 

-168:^5(2, -1) - 2 JP-5(3, -3) + 180:^5(3, -2) (A.8) 

, , 1 1433 1 , , 461 , , 85 , 

(A.lOj 

For positive q we have 

427 11 5 

19 , V 1 , V 13 ^ V 35 , V 85 

^^s{l, 0) - 7^^5(1, -1) + 7^^5(2, 0) + —^5(2, -1) + -7r7^5(2, -2) 



768 ' ' 9216 ' ' 48 ' ' 288 ' ' 2304 

^^(3, -4) - ^^5(3, -3) - ^^5(3, -2) (A.ll) 



6144 ' ' 512 ' ' 2304 

= 96^ " 1024^ " 1105 9 2 007r2 ~ 30-^'^°' + 2880-^'^°' 
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_i) + ^^,(2, 0) - Jl^^,(2, -1) - ^ 

276480 ' ^ 11520 ' ^ 276480 ' ^ 13824 

4433 _ , 437 _ 161 _ 

- J^A 3, -2 + -^^JT. 3, -3 + -^.r^TAZ, -4 



(A.12) 



4433 _ , 437 _ 161 _ 

^ . X 1 3 214253 13 ^77^ 

" ~ 1024^ ~ 58982400^ ~ 320 + 3840 

47 149 391 169 

+ ^^5(0, 1) + —J^sih 2) + ^J^sil, 1) + — ^J^5(l, 0) 

92160 ^ ' ^ 960 ^ ' ^ 7680 ^ ' ^ 294912 ' ' 

7 37 323 119 

6139 399 49 
^!^jr (3 _2) + Jl^jr (3 _3) + _^jr (3 _4) (A.13) 

491520 ^ ' ' 40960 ^ ' ^ 32768 ' ' ^ ' 

1 1 43 230795603 

' ^ ~ 1927r2m6 ^ 38407rW + 691207r%i2 ~ 2774532096007r2 

22817 ^ 3529 ^ 509 ^ 533 ^ 

"1935360-^^^°' ~ 3870720-^^^°' + 1720320-^^^°' ~ 6A5U0^'^^^ ^> 

164753 _ , , 2823307 _ , , 19 _ , , 25471 _ , 

+ ^51, 1) 1, -1 J^5(2, 

5160960 ^ ^ 990904320 ^ ^ 55050240 ^ ^ 860160 ^ ^ 

8917219 323 17909419 

+ 1486356480-^'^ ' + 2752512 ' ' ~ 2972712960-^'^ ' 

1083_ _ ^JF5(3, -4) (A.14) 

9175040 ^ ' ' 1048576 ' ^ ^ ^ 

r 1 19 1 " 

^,(0,4) = (1 - 51ogm^) [^g^ - ^g^g^ + g^Ye^ 

31 13 1 

-^^.(0,3)-— ^,(0.2) + — ^.(0,1) 



+s 



t4 ^ ' 1152 ^ ' 9216 ^ ' 

5 61 80989 347 

u 5767r2m4 + 184327r2m2 + 884736007r2 + 1440 ' 

_^^,(0, 2) + ^^^,(0, 1) - 2) + ^^^^(1. 1) 

2560 ^ ' ' 184320 ^ ' ' 2880 ^ ' ' 23040 ^ ' ^ 

415 , 23 _ , , 329 _ 13283 _ 

-147456^'^^' °) + 4423680^'^^' " TT520^'(2' + Tl05920^'(2' 

"22Tl84-^'^^' "^^ " 22TT840-^'^^' "^^ + 245760-^'^^' "^^ + 589824-^'^^' 

(A.. 

2n r 1 1 415 11 " 

J^5(0,5) - (1-Mogm ) [ 1927^2^6 + issGTr^m^ ~ 44236871^ + 4423687r2. 

+ J^^,(0, 3) - ^^^5(0, 2) + ^^^,(0, 1) 
13824 ' ^ 110592 ' ^ 884736 ^ ' ^ 

r_ 7 _ 101 1207 6492275537 

+ [~23 047r2m6 ~ 92 1607r2m4 + 176947271^ + 99883155456007r2 
634603 2623 61477 172157 

7741440-^'^°' " 258048-^'^°' + 371589120 •^'^°' ^> ~ 2580480-^'^^' 
515491 778597 7297 

61931520 ' ' ~ 1321205760 ' ~ 3567255552 ' ~^ 
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+ 
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57917 32033909 _^ 620245 

"10321920 ' ^ ^ 17836277760 '~ ^ ^ 891813888 ' ~ ^ 
_^05348^ _ 138643 _ 36485 1 

7134511104 '^^ ' ' 198180864 ' ' 339738624 '^^ ' ^ ' 

(0 6) - (1 5\o m^) \ ^ 1 1 6209 

I , j - ( ogm ) [3207^2^8 + 38 407r2m6 ^ 20 4 807r2m4 353894407r2m2 

1091 1 1429 4373 ^ 401 

2831155207r2j 221184 ' ^ 8847360 ^ ' ^ 70778880 ^ ' ^ 

9 77 709 24257 

.~64007r2m8 ~ 2304007r%^ ~ 860160071^ + 1189085184^1^ 
903503350441 74410051 1058653 

47943914618880007r2 ^ 3715891200 ' ^ ~ 464486400 ' 
982249 16644221 ^ 2348599 

59454259200 ' ^ 1238630400 ' ' 9909043200 ' ' 

791158043 3838603 , 3179779 , 

7-^ f 1 0^ 7^ (^ 1 1 -I 7^ (Q. ol 

5707608883200 ' ^ 8561413324800 ''^ ' ' 4954521600 ' ' 

466701533 (2 _i) 65256251 ^2 _2) 9795127151 _2) 
"1223059046400 '~ ^^ 428070666240 '~ ^ ^ 17122826649600 ' ~ ^ 

72933457 ^ /o _on _ 3838603 (o _a{\ 
"475634073600-^'^^' "^^ " 163074539520-^'^^' " ^ " 
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